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Abstract. We discuss the possibility that we could obtain some hints of the heavy physics 
during inflation by analyzing local features of the primordial bispectrum. A heavy scalar 
field could leave large signatures in the primordial spectra through the parametric resonance 
between its background oscillation and the fluctuations. Since the duration of the heavy- 
mode oscillations is finite, the effect of the resonance is localized in momentum space. In this 
paper, we show that the bispectrum is amplified when such a resonance occurs, and that the 
peak amplitude of the feature can be O(10 1-2 ), or as large as O(10 5 ) depending on the type 
of interactions. In particular, the resonance can give large contributions in finitely squeezed 
configurations, while the bispectrum cannot be large in the exact squeezed limit. We also 
find that there is a relation between the scales at which the features appear in the bispectrum 
and the power spectrum, and that the feature in the bispectrum can be much larger than 
that in the power spectrum. If correlated features are observed at characteristic scales in the 
primordial spectra, it will indicate the presence of heavy degrees of freedom. By analyzing 
these features, we may be able to obtain some information on the physics behind inflation. 
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1 Introduction 

Primordial non-Gaussianity is a powerful probe to discriminate inflationary models [1]. In the 
simplest setup, i.e. single-field slow-roll inflation models with a canonical kinetic term and 
Bunch-Davis initial conditions, the primordial fluctuations are predicted to be approximately 
Gaussian-distributed [2]. More generally, for the models of so-called single-clock inflation, 
where there is only one relevant degree of freedom, it is known that there is a consistency 
relation relating the bispectrum in the squeezed limit to the tilt of the power spectrum [3-8] . 
Thus, a non-Gaussian signal would enable us to narrow down the possible models of inflation. 

To predict the expected non-Gaussian signal from a given model correctly, we need 
accurate knowledge of the inflationary dynamics as well as the Lagrangian. For example, even 
in the "single-field" inflation models, where there is a single light scalar field, the consistency 
relation is not always satisfied. It has been argued in Ref. [9] that when the background 
evolution is provided by a non-attractor solution, the consistency relation is violated while 
maintaining a scale-invariant power spectrum. To probe the model of inflation, we need to 
understand well what can happen under a given model. 

In general, there could be many scalar degrees of freedom other than the inflaton, even 
in single-field inflation models. In a model embedded in supergravity or string theory, for 
example, such degrees of freedom may appear as moduli fields, Kaluza-Klein modes, or the 
scalar super symmetric partner of inflaton. Usually, the scalar fields are very heavy, m S> H, 
so that they are assumed to be stuck in their minima, and the model is treated effectively as 
a single-field model [10]. Recently, however, it has been pointed out that a heavy scalar field 
is not necessarily stuck in its potential minimum during inflation [11-20]. It can be displaced 
from its minimum due to the centrifugal force generated by a turn in the inflaton trajectory. 
When the turn is very sharp, even oscillations in the heavy direction can be excited [19, 20]. 
Other possibilities are that oscillations can be excited when the heavy scalar field becomes 
momentarily light /tachyonic during inflation or that they can be excited at the beginning 
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of inflation [21], which is natural in the case that inflation occurs after tunneling from a 
neighboring minimum [22-27], for example. 

In the previous paper [28], we discussed the possibility that excited oscillation of heavy 
modes can leave non-negligible signatures in the power spectrum through derivative cou- 
plings, without spoiling inflation (see also Refs. [29] for models with only gravitational 
couplings). We saw that the primordial fluctuations can be enhanced deep in the horizon, 
k/a ~ m S> H, by the parametric resonance with the excited oscillations. In this paper, 
we estimate the resonant feature in the bispectrum within the same setup. This gives an 
example where resonant non-Gaussianity [30-34] is produced in a realistic model. Since the 
duration of heavy-mode oscillations is finite, the effect of the resonance is localized in mo- 
mentum space. We will show that a large feature could be induced in the bispectrum even 
when the feature in the power spectrum is too small to be observed. We also investigate the 
behavior of the bispectrum in the squeezed limit. As in the case that the consistency relation 
is satisfied, it can be shown that the bispectrum in the squeezed limit cannot be large unless 
the modification to the power spectrum is large. However, we will also show that, for finitely 
squeezed configurations, the bispectrum can be greatly enhanced in comparison to the usual 
single- field inflation model (a model with a single light scalar field). Since we cannot ob- 
serve the exact squeezed limit in actual observations, the enhancement could be practically 
important, as recently pointed out in Ref. [35], where the resonant non-Gaussianity from 
a modulated potential is considered as a concrete example. The detection of such features 
would therefore indicates the presence of heavy degrees of freedom during inflation, and by 
analyzing them we may hope to improve our understanding of the physics behind inflation. 

The organization of this paper is as follows. In §2, we briefly review the model pre- 
sented in Ref. [28] to realize an efficient enhancement of the fluctuations in the inflaton field 
and estimate the feature induced in the power spectrum. In §3, we discuss the resonant 
enhancement of the bispectrum. Finally, we provide a summary of this paper in §4. 



2 Overview of the Model 
2.1 The Model 

First, we briefly review the model presented in Ref. [28]. We consider a model with a heavy 
scalar field with mass m S> H, x(= ^ )i which derivatively couples to the inflaton field, 

<j>{= <^> (1) ), as 



S m = - f dx 4 ^P(X IJ , 4> K ), (I, J, K = 1, 2) 



(2.1) 



dx \f—g 



\(d4>f + V(cf>) + \{d X ) 2 + T ^x 2 + K n + K d + {X(b 2 /AlX 3 /A 8 d ) 



(2.2) 



with X IJ = —d fl 4 >I d tJ- (t) J /2. Here, we have assumed that the derivative couplings at the 
leading order are given by, 

A, 



K„ 



2A, 



and 



K d 



(2.3) 
(2.4) 
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with the dimensional parameters A n ,A^ and the dimensionless parameters A n ,A^j = 1)2). 
The derivative couplings in Eq. (2.2) generally appear in the action from the effective-field- 
thoery point of view [36, 37] . The derivative couplings K n and provide the most general 
couplings between the inflaton field and the heavy scalar field at the leading order in 1/A ?t 
and in a model with the parity symmetry, <j) — > —(f), and the shift symmetry, (j) — > <j) + c. 
The approximate shift symmetry is usually assumed to ensure the flatness of the inflaton 
potential. In addition, we have assumed the parity symmetry to forbid the kinetic mixing, 
which obscures the difference between the light and heavy fields. Thus, we first focus on the 
couplings K n and Kj, though brief comments are presented in §3.4 on expected features in 
the bispectrum when a symmetry-violating interaction becomes relevant. In general, higher- 
order terms in A n and A^ are also expected to appear. To ensure that contributions from 
these terms can be treated perturbatively, we assume here that the background fields satisfy 
the following conditions, 

X < A n , 0, X < A^. (2.5) 

We have also suppressed the terms (dcfi) and (dx) 4 in Eq. (2.2), because they have a little 
effect on our analysis for the background evolution and the power spectrum. The higher- 
order terms and the self interaction (dcj)) 4 can be important for the bispectrum. Features in 
the bispectrum induced by these interactions are discussed in §3.4. We can also find many 
concrete models which have the couplings K n and Kj, such as DBI inflation [38]. In that 
case, the conditions (2.5) are not mandatory. 

Since x can generally decay into other particles, we assume that x decays with a rate 
r, which satisfies H <C T <C m. To ensure the slow-roll inflation, we further make two 
assumptions. First, the inflaton potential is sufficiently flat, 



where 



ey <C 1, |wl < 1, (2-6) 



Ml fv'Y 9 V" , . 

zv=^{v) ' W^-, (2.7) 



are the slow-roll parameters. Here, M p = 2.4 x 10 18 GeV is the reduced Planck mass. Second, 
the heavy scalar field x is subdominant, 

fx « (2-8) 

where 

_ Px _ X 2 + m 2 X 2 
Ix - p - 6M|iJ 2 ' 1 ' 

is the fraction of its energy density to the total one. 

Provided that the conditions (2.5), (2.6), and (2.8) are satisfied, the background evolu- 
tion of the inflation field is given by the slow-roll solution, 

^(t) * ~, (2.10) 
where tt^ is the conjugate momentum for the inflaton field, 

7T = Z 2 J, (2.11) 
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Figure 1. The fluctuations in the inflaton field can be amplified through the parametric resonance 
with the background oscillation of the heavy scalar field (2.13). As depicted by the thick line, a mode 
is redshifted by the cosmic expansion. Since only the modes that have crossed the resonance band 
during the oscillation are amplified, the features in the primordial spectra are localized in momentum 
space. 



with 



l + A * + (A, 1 + A d2 )^ 



(2.12) 



and = 2? /a 3 . The heavy scalar field oscillates with a frequency of the mass scale m, 



X(t) ^ Xoe rt cos{mt)9(t), 



(2.13) 



with 9{t) being the Heaviside function, where we have assumed that the oscillation is excited 
instantaneously at t = 0. The fluctuations are enhanced through a resonance with the 
background oscillation of the heavy scalar field (2.13). We showed the basic picture of the 
resonant enhancement in Fig. 1. 

For a later convenience, we introduce the parameters which represent the magnitude of 
the derivative couplings, 

Xo 



A, 



Qdl = Arfi 



2 2 

m Xo 
2A^ ' 



Qd2 — Ad2 



2 2 

Xo 
2A^ : 



(2.14) 
(2.15) 



which are related to the g-parameters introduced in Ref. [28]. We will also use q to represent 
the order of the (/-parameters above. Note that the g-parameters are less than unity under 
the conditions (2.5). They can be expressed in terms of the energy fraction of the heavy 
scalar field, f x , as 



V2\ r 



inf 



A, 



f? 



for the couplings K n and 



Qdl — A<fl fx,t=o, Qd2 — Ad2 



X,*=0' 



Eini 



An 



fx, 



t=0, 



(2.16) 



(2.17) 



for the couplings , where E- m { is the energy scale of the inflation, E- m { = (SMrH ) 4 — V * . 
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2.2 Power Spectrum 

Before investigating the features in the bispectrum, we first estimate the effect of the reso- 
nance on the power spectrum. To obtain an analytic expression, here, we use the perturbative 
method assuming that the g-parameters are sufficiently small. 

To eliminate the gauge degrees of freedom in the metric fluctuations, we use the flat 
gauge, where the spatial metric becomes a 2 5ij. In the flat gauge, we can make some simpli- 
fications in estimating the power spectrum. First, the metric fluctuations can be neglected 
because the resonance takes place deep in the horizon. Moreover, we can also neglect the fluc- 
tuations in the heavy scalar field, 5\- This is because 5\ oscillates with a frequency ui > m 
and then does not satisfy the resonance condition. In addition, the mass term suppresses the 
amplitude of Sx at k/a ~ m compared to that of the inflaton field. 

Neglecting the fluctuations in the metric and the heavy scalar field, the second-order 
action for the fluctuations in the inflaton field, <p, can be written as, 



S 2 ~ J did 3 * | 



c 2 s (V<pf/a 2 



(2.18) 



where 



1 



2A^ 



+ o 



x_ 



(2.19) 



is the speed of sound for the fluctuations in the inflaton field. Note that we can neglect terms 
like <p 2 /A d in the action (2.18), which could be induced by the self interaction (9^) 4 /A^, 
because they are higher order in the (/-parameters. 

Introducing the canonically normalized variable v = z^tp, the action (2.18) can be 
rewritten as, 



So 



dtd d x 



and then the full Hamiltonian as, 



v 2 



v 2 + 



dV 2 



c 2 V 



2* 



Z<f, 



(2.20) 



(2.21) 



Substituting the background solution (2.13), the quantities c 2 and z^/z^ can be estimated 
as, 



l + ^V 2r *cos(2mt), 



m 



q n e rt cos(?7ii) + (q dl + q d2 )e 2Vt cos(2mt) + O ( — J 

V ml 



(2.22) 
(2.23) 



In the following analysis, we neglect the 0{H/m) terms in Eq. (2.23) assuming that q is 
larger than H/m. Taking the free Hamiltonian as the terms independent of the derivative 
couplings in Eq. (2.21), the interaction Hamiltonian is given by, 1 



1 Note that infinitely many terms would appear in the interaction Hamiltonian if we define the free Hamilto- 
nian as terms independent of the derivative couplings without introducing the canonically normalized variable 
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H 



(2) 



m 



2(2vr) 3 

where we have moved to Fourier space, 

1 



d 3 k crV<;-k, 



v(x) 



(2nf 



d k v^e 



ik-x 



(2) 

Here, the coefficient Ci is given in terms of the g-parameters as, 



c (2) ^ C (2;n) e - ri CQs(mt) + C (2; d ) e -2K cog(2m ^ 



where 



(2;n) 



1 



k 



1 



Qdl- 



(2.24) 
(2.25) 

(2.26) 

(2.27) 
(2.28) 



2 \am / 

Hence, using the in-in formaiism, the correction to the power spectrum can be estimated as, 



at the leading order in the g-parameters, where the fields in LHS are operators in the Heisen- 
berg picture while, in RHS, they represent operators in the interaction picture. For brevity, 
we denote them in the same way since it will be clear which they represent from the context. 
In terms of the dimensionless power spectrum, V<p, defined through, 



^k(%k'W> = (2vr)^(k + k / ) 



a2tt 2 7? 



k 3 



the correction is expressed as, 



r(2) 



where 



(2) 



-2m 2 Re 



u k (t 



>\2 



(2.30) 



(2.31) 



(2.32) 



Here, the estimation time t should be taken sufficiently after the decay of the oscillation. 
The function u k is the mode function for v^, which is defined through, 



^k(t) = u k (t)a k + u* k {t)a. 



t 



(2.33) 



for the creation/annihilation operators. Though the state could be excited depending on the 
excitation mechanism of the oscillation, we simply assume that the state was in the vacuum 
state at the excitation of the oscillation. Even if the state was excited, the resonance is not 
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spoiled except for some specific excited states. Assuming that the state is not excited, the 
mode function is provided by, 




«*(*)= \hr [ l -TZ (2-34) 



where r is the conformal time. Because the resonance occurs deep in the horizon, we can 
approximate the mode functions in the integral (2.32) as, 

u k {t) ~ J^e~ ikT - (2-35) 

Hence, 1^. can be approximated as, 

4 2) - " r dz> (—) C k ] sin( 2 kr'), {z = mt). (2.36) 
Jo \amj 

The resonance occurs between the factor sin(2/cr') and the oscillatory components, cos(mi') 

(2) 

andcos(2mt'), in ' (see Eq. (2.26)). If the negative energy mode exists, it would introduce 
an additional oscillatory component in Eq. (2.36) with a constant phase shift. However, this 
additional component does not spoil the resonance unless the amplitude of the negative 
energy mode and the phase shift have specific values for the scales of the resonance. 

(2) 

Here, we evaluate /£ for the couplings K n . By replacing the parameters (q n , m/2, T/2) 
by (<7d, m, r), we can obtain a similar result for the couplings K^- Using the product-to-sum 
identities for the trigonometric functions, we have two oscillatory functions in Eq. (2.36). 
The integration has a resonant contribution from one of the two in the interval around the 
stationary point of the phase function, #( 2;n ) = mt — 2kr. Taking the derivative of the phase 
function 6^ 2 ' n \ we can find that the resonance occurs at time £* where, 

km , 

(2.37) 



a(t*) 2 ' 

Hence, the resonance occurs AiV = ln(m/2H) e- folds before the horizon crossing of the 
modes. Solving the above equation, the value of z at the resonance is estimated to be 



m 



2k \ 



z* ~ - In , (2.38) 

H \a^m J 

where ao is the scale factor at the onset of the oscillation. Then, the integral can be evaluated 

as, 2 

if cC^Im/W, (2.39) 

where 

jW= f Z dz'e-^ z ' +i ^'~ 2kT '\ (2.40) 



2 Note that we have used 2k for the argument of for later convenience. The factor 2k in Eq. (2.40) 
corresponds to k + k' if we don't take into account the delta function in Eq. (2.30). Similarly, the scaling 
function 1^ with K = k\ + &2 + &3 appears for the bispectrum. 
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and c£ 2 ^ is the coefficient function (2.27) at the resonance, t = t*. The duration of the 
resonance can be roughly estimated by Az = m/\/~0, which is estimated to be 

Az^ = \/^/H, (2.41) 

for the couplings K n . Approximating the oscillating factor e l ( z '~ 2kr ') [ n Eq. (2.40) by a 
top-hat function, then, we can roughly estimate the contribution to I^u from the resonance 
through the couplings K n as, 



An) n (2;n) . J2;n) [™*<f>,*.+V%) , _ r 



max(0,2»-Jf) 



^(2;n) . J2;n) An) 
°jfc,* bmf7 * J 2fc,approx 



2 ^ | H P %/Tlr, 



■ /i(2;n) 



2m 2k \ H 



r v ao m 



sinh 



' Hm 



8111 



,(2;n). 



In ( ] < 

1 aom ' 



K < ln ( _2*_ , < 
m \ aom ' 



ln(-^-) > /H 

1 aom 1 



(2.42) 
(2.43) 

(2.44) 



where oi 2 '^ is the phase function at the resonance, 



i(2;n) 



77i 

7F 



In 



2/c \ 



aom 



+ i 



(2.45) 



Hence, the power spectrum has a peak at k p /ao ~ me s J u l m /2 with an amplitude, 



A^)(fc p )/P, 



m 



1 — e 



(2.46) 



In particular, when the decay rate is small, T <C \/Hm, the peak amplitude becomes, 

AV^{k p )/V^2q n ^ (2.47) 

which is of the order of 2q n Az^ 2 ' n \ 

In Fig. 2, we showed the scaling function (Eq. (2.40)). The function nas a peak 
at kja^ra ~ 1. As can be seen in the figure, the approximate function ^™ a pp rox introduced in 

Eq. (2.44) gives a good approximation of . 

Once the correlation function for the inflaton field are obtained, we can calculate those 
for the comoving curvature fluctuations £ by using the gauge transformation. 3 After the 
oscillation damp out, the gauge transformation is given as usual, 



r H 
C = ~rf, 



(2.48) 



3 Note that the notation 1Z is also used to denote the comoving curvature fluctuations. The variable £ here 
corresponds to 1Z in Ref. [39], for example. 
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Figure 2. The scaling function IV*' (Eq. (2.40)), which determines the scale dependence of the 
resonant enhancement. The thick line represents the approximate function, If. a pprox - The function 
has a peak at k/ (aom) ~ eV ff / m . In the plot, we have set T = 5-ff and m = 10 3 if. 



at the linear order. At the non-linear level, in general, we can obtain the correlation functions 
of ( by using 5N- formalism [40-43] from those of <p on superhorizon scales. Since £ is 
proportional to tp at the linear order, the correction to the power spectrum for £, AV^/V^, 
is provided also by Eq. (2.31): 

AVdk)/V c = /f * cgj»> S in#"^ approx . (2.49) 
Hence, the peak amplitude is given by Eq. (2.46) (or Eq. (2.47)): 

AV^\k p )/V c ~ ^ (l - eTviw) g n (2.50) 
~ 2Qn \flj (for T< Vim). (2.51) 

We close this section with a comment on the perturbativity. Even when the ^-parameters 
are smaller than unity, the higher-order terms in the expansion with respect to the interac- 
tions can be comparable to the leading term when the duration of the resonance At is 
sufficiently long. This is because the resonance coherently accumulates the effects of the 

(2) 

interaction. Then, the order parameter of the expansion is given by Hj At, which could be 
large even if q is very small. The perturbative method is valid only when the correction to 
the power spectrum (2.47) is less than unity. When this condition is not satisfied, we should 
solve the full equation of motion numerically to get the mode function as done in Ref . [28] . 

Other than the resonance, the derivative couplings could induce the loop corrections to 
the power spectrum and the Green function. To ensure that the corrections are suppressed, 
the model (2.2) should have a cutoff at the energy scale A cut = min(2-7rA n , (2-k) 1 / 2 kd). 4 
Hence, to perform the calculation perturbatively at the resonance scale, k/a ~ m, the mass 
scale should not exceed the cutoff scale, m < A cut . 

4 Here, we have assumed that the coupling constants of the derivative couplings, A, are of the order of unity. 
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3 Features in the Bispectrum 



Next, we investigate the effect of the resonance on the bispectrum, 

(C kl (t)Ck 2 (t)Ck3(*)) = (2vr) 7 5 3 (k 1 + k 2 + k 3 )P c 2 -§-2 • (3-1) 

K l K 2 K 3 

After seeing which cubic interactions are important for the resonance, we estimate the cor- 
rection to the bispectrum perturbatively by using the in-in formalism. 

We use the dimensionless quantity B^ to represent the amplitude of the bispectrum, 
which is made dimensionless by the factor k\k\k\ as in Ref. [44]. Note that this normalization 
factor differs from that in the definition of the local- type /nl- In addition, we use the 
uncorrected power spectrum to define the amplitude Bq. 

3.1 Interactions Relevant to the Resonance 

Here, we discuss which cubic interactions play an important role in the resonance. First, 
they should obviously contain oscillatory components. Secondly, the interactions with more 
derivatives become important for the resonance. This is because the resonance occurs deep 
in the horizon, k/a ~ m 3> H , where the fluctuations oscillate with high frequencies. Note 
that the fluctuations in the heavy scalar field can be neglected at the leading order as in 
the case of the power spectrum. Hence, we do not need to consider the mixing between the 
inflaton field and the heavy scalar field. 

Taking into account the above points, we can pick up the relevant interactions as (see 
Appendix A), 



„(3) _ a 3 









r %\ 


j d 3 x < 






l 2 . 



2 



(3 - 2c 2 )0 2 + |(V^) 2 

+ [z\ + 2(q dl + q d2 )e- 2Tt sin 2 (mi)] V^V^ViW j. (3.2) 

Note that we need to include the gravitational couplings to obtain the interactions (3.2) 
because the couplings K n and K d do not contain cubic interactions such as [d\ • dcft)(d(j)) 2 
due to the parity symmetry. This is why the interactions (3.2) are suppressed by the Planck 
scale, M p . However, direct cubic interactions exist if the self interaction (9^) 4 or higher-order 
interactions are introduced, though we have neglected them in the previous section because 
they have a little effect on the background evolution and the power spectrum. They are also 
expected to arise if we relax the requirement of the parity symmetry. We will discuss these 
possibilities in §3.4. 

In terms of the canonically normalized variable v = z^tp, the interaction Hamiltonian 
can be written as, 



Hf ] = ^—^ J^. [ d^fedVCki +k 2 +ks) 

a2Mj27r) e V 2 J 



2^(30) 



a2M p (27r) £ 

+ mC7 k!k 2 k3^ k i K 2 v k3 y + C^Wk^ik,] • (3.3) 
Here, each coefficients are decomposed into the K n - and i^-coupling parts as, 

cS 2 k3 - Cgia^ cos(mt) + Cg^e- 2 ™ cos(2mt), (i = O, II), (3.4) 
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and 



where 



c£L - cg^e"" Mmt) + Cg^e" 2 " sin(2^), (3.5) 



^,(30;?!) gn ■ k 2 , . 

k lk2 k 3 ~ 4 (flm) 2 > ^ 



and 



kik 2 k 3 2 y2 k\ 

r (3ii;n) _q n (1 _ ki • k 2 , 
k!k 2 k3 - 2 1^2 k\ 



r {30;d) _ Qdi + 3g<a ki • k 2 . . 

k lk2 k3- 8 (am)2 > ^- y i 

^(3/;d) _ qdi + qd2 ki-k 2 \ , , 

C k lk2 k 3 = 2 U " "fcf - J ' { ] 



n (UI;d) _ Idi (1 . 3ki • k 2 \ 3g d2 / 1 ki ■ k 2 . 
C k lk2 k 3 = ^2 + _ fcf - ) + ~\2~ ' • 



Here, we have extracted only oscillatory components from the interactions (3.2) and neglected 
the higher-order terms in the g-parameters and H/m. 

3.2 Resonant Enhancement of the Bispectrum 

Provided the interaction Hamiltonian (3.3), we can estimate the correction to the bispectrum 
by using the in-in formalism, 

W*WtWt)> - i fdt' <[tfiV)^kJ%k 2 (%k 3 (0]>, (3-12) 



at the leading order in the ^-parameters. To see the contributions from the resonance, it is 
sufficient to consider the linear gauge transformation (2.48), £ = Hip/<f>. Then, in terms of 
the dimensionless bispectrum (3.1), the correction can be estimated as, 



^C-f'gU,. (3-13) 



Here, 



/ k 1 ) k2 k 3 - ** 



Z dz' a 



H In klk2k3 



-iKr' 



+ (5 perms), (K = k x + k 2 + k 3 ), (3.14) 



where the coefficient function C^V ,, is defined as, 

KiK 2 K3 ' 

r (3) = r (30) _ . / k 2 + k 3 \ (37 ) _ / k2_\ f k3_\ (3J7 ) 

°kik 2 k3 - °kik 2 k 3 M am I u kik 2 k 3 I am J \ am J U kik 2 k 3 - ^- i0 J 

We have also introduced the averaged slow-roll parameter ey by, 

ey = ^£v,*ev,c (3.16) 
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where ey,* and ey )C are evaluated at the resonance and the horizon crossing, respectively. 
In deriving Eq. (3.14), we have used the subhorizon approximation for all mode functions, 
u ki{t') (* = I? 2, 3), assuming that all modes were in the subhorizon regime at the resonance, 
ki/a^H > 1. This approximation is not valid when we consider the squeezed limit, where a 
single mode can be outside of the horizon when the oscillation is excited. We consider the 
squeezed limit separately in the last part of this subsection. 

Since the integral has a similar form as Eq. (2.36), we can estimate it as in the previous 
section. The oscillatory components, e imt ' and e 2%mt \ in C^^, k (see Eqs. (3.4) and (3.5)) 
resonate with the factor from the mode functions, e~ lKr ' . The resonance occurs in a similar 
way as in the power spectrum, where the oscillatory components in the interactions resonate 
with the factor from the mode functions e~ 2lkr ' . Here, we evaluate ^jj- je k ^ or ^ ne couplings 
K n . By replacing the parameters (q n ,m/2,T /2) by (qd,m,T), we can obtain a similar result 
for the couplings K^. The phase function #( 3;n ) = mt — Kt has a stationary point at time i* 
where 



K 

a(U) 



m. 



(3.17) 



Then, using the function (2.40) with K in the argument, the integral can be evaluated as, 



(3) 



kik 2 k 3 



~ Re/ 



K 



a 



(30) _ k 2 + k 3 ( 3 j) 

kik 2 k 3 ,* j£ Ly kik 2 k 3 ,* 



K [k k lk2 k 3 ,* 



m 

Jj, (3-18) 



where the coefficient functions are evaluated at the resonance, t = t*. The duration of the 
resonance is given by Eq. (2.41). Note that the bispectrum has contributions from the 
subhorizon regime [30-34] in contrast to the cases without a resonance. Then, 1^ ^ can 
be approximated to be, 



T (3\n) /i(3;n) T (n) 

i, , , ~ COS Ui I\r 
kik 2 K3 * K,approx 



a 



(SO) 



kik 2 k 3 ,* 



k2 + h n {31) 



K 



kik 2 k 3 ,* 



h\ c (3ii) 

J£ J kik 2 k 3l * 



m 
H 



+ (5 perms), 

(3.19) 



and then the correction to the bispectrum (3.13) can be roughly estimated as, 



~ 4 yjjJQnCOSV* J ^ ia pp rox ,5 ki,k 2 ,k 3 , 



(3.20) 



where 

<Ski,k 2 ,k 3 = Q n 



r (30) _ fc 2 + k 3 (37) _ (tl) (tl) C {3H) 

kik 2 k 3 ,* kik 2 k 3 ,* \ K ) \ K I k ik 2 k 3 ,* 



+ (5 perms), 

(3.21) 



for the couplings K n . Therefore, the amplitude of the correction to the bispectrum can be 
written in terms of that to the power spectrum (2.50) as, 



(n) ^f™)^( 

t 4 \HJ Vc 



(») 



"<5ki,k 2 , 



k 3 - 



(3.22) 
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1.0 



Figure 3. The resonant enhancement of the bispectrum through the couplings K n . In the left panel, 
we showed the shape function (3.21). Along the solid line, the bispectrum scales as shown in the right 
panel. The lines indicate the same in Fig. 2. The bispectrum has large values along the dashed line, 
K = k\ + &2 + /C3 = 2fc p , where k p is the peak scale of the feature in the power spectrum. Here, we 
have set r = 5H, m = 10 3 H, q n = 0.1, and ey = 0.01. If other interactions exist, larger enhancements 
can be realized (§3.4). 



As is clear from the expression (3.22), the correction to the bispectrum can be much larger 
than that to the power spectrum thanks to the large factor m/H. In Fig. 3, we showed the 
shape and scaling of the correction to the bispectrum. The bispectrum has large values for 
K = 2k p , where k p is the peak scale of the feature in the power spectrum, k p /aom ~ 
Therefore, the characteristic scale of the feature in the bispectrum is correlated with that in 
the power spectrum. 

To see the magnitude of the amplification, here, we compare the result (3.20) with the 
equilateral-form and local-form bispectra at the specific configurations. 



- Equilateral configurations {k\ ~ /c 2 ~ /c 3 ~ 2k p /3) 
The amplitude of the correction is estimated to be 



^.equilateral - IQ \ H J p ( 

- lT qn (9) § (for r <<; ^^ } ' (3 ' 24) 

at the equilateral point k\ ~ k2 — k$ ~ 2k p /3, while for the equilateral-type /nl, 

^Ci e q uuat cral = Yq-^nl ■ (3.25) 



(n) 



- Squeezed configurations (k^ <C k\ — k<i ~ k p ) 

In the finitely squeezed configurations at the specific scale, ks k\ — ki ~ k p but 
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&3 > a *H ~ (H/m)K, we obtain a non-vanishing contribution, 



AB 



(n) _ 7ey /m\ AP, 

(,not so squeezed -^g \JJ / 'P 



(n) 



while for the local- type /nl, 



R ^ # l oca l ( k p 



(3.26) 
(3.27) 

(3.28) 



In terms of the bispectrum, 

(2vr) 3 5 3 (\>) P c (fc x )AP< 



7ey /m 

IT If 



(3.29) 



where P<^ is the power spectrum for £, P^ = 27r 2 P^/A; 3 . On the other hand, in the standard 
single- field slow-roll inflation models, the bispectrum for the squeezed configurations is of the 
order of the slow-roll parameters: 



(CidWCk^RksW}^ 



(2vr) 3 <5 3 QTiq P^P^) 



(3.30) 



where n s is the spectral index of the power spectrum. Hence, the excitation of a heavy scalar 
field can give a non-negligible effect on the bispectrum for the configurations, (2H/m)k p <C 
&3 <C fei — — fcpj as recently pointed out in Ref. [35]. If the peak scale k p is not so large, 
the scales k <C (2H/m)k p could be in the superhorizon regime at the present time. In that 
case, these contributions could provide practically the most important one in observations 
of the bispectrum in the squeezed limit. Note that the contribution (3.29) corresponds to 
the 0(ks/kp) 2 correction in Eq. (3.30). This correction can be large in the resonant cases 
because it appears in the combination (k^T*) 2 ~ (k^/2kp) 2 (m/H) 2 , which is not necessarily 
small even for the squeezed configurations k% <C k\ ~ ki- 

As mentioned in the text below Eq. (3.14), our subhorizon approximation is not valid 
in the squeezed limit where k^ <C (2H /m)k p , k\ ~ k<i ~ k p . In this case, the mode function 
Uk 3 and its time derivative iik 3 should be replaced by those with the extra factors i/kr' and 
l/(fer') 2 , respectively. Hence, the contributions to the squeezed limit can be estimated as, 



AB 



(«) 

£ squeezed 



e v q r 











f kp 


\ fm 







and the bispectrum behaves as, 

(Ck 1 (<)Ck a (*)Ck3(t))*3<fc 1 ^=*p ~ (2^) 3 ^ 3 (E k j°c(*i) A ^ (n) (*3)^- 



(3.31) 
(3.32) 

(3.33) 
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Therefore, the bispectrum scales as in the usual case (Eq. (3.30)) and cannot become large 
unless the modification in the power spectrum is large. 

As discussed in the previous section, we need to require the conditions AV^ /V^ < 
1 and m < 2irA n for the perturbativity. The correction has the maximum value when 
these conditions are saturated. In this case, both energy scales m and 2nA n are given by 
0(10 3 )(/ Xi j=o/ey) 1 ^ 3 H . Then, our perturbative calculation is reliable if 

AflW ~ e v Q < 0(1O 3 ) x (4fx,t=o)^. (3.34) 

Here, we have used the observed value for the amplitude of the power spectrum, = 
0(1O -9 ) [45]. Similarly, in the case of the couplings , we can find the upper limit for the 
perturbativity as, 

Atff ~ e v (g) < 0(1O 3 ) x (4&m>)*. (3-35) 

The energy scales m and (27r) 1//2 Ad are given by O(10 3 )(f x>t=0 /ev) 2/7 H when the inequality 
is saturated. 

Therefore, the correction can be ©(lO 1-2 ), which is much larger than that expected in 
the standard single-field slow-roll inflation models, even when the correction to the power 
spectrum is small. More larger non-Gaussianity can be induced when the correction to the 
power spectrum is measurably large. 

3.3 Perturbativity 

As in the case of the power spectrum, the higher-order terms could be comparable to the 
leading term even when the q-parameters are small. Here, we discuss the validity of the 
perturbation for the cubic interactions (3.2). 

First, we estimate the magnitude of the interaction Hamiltonian (3.2) in the resonance 
regime from the dimensional argument. The interaction Hamiltonian (3.2) can be estimated 
as, 

„(3) aq^/e^ 8 3 

where we have assumed that the wavenumbers have the same order k. The fluctuations ^ 
can be estimated as, 

^~MW- (3 - 37) 

Here, we have replaced the creation/annihilation operators by {2n/k)^, taking into account 
their commutation relation [ak) a k' ] = (27r) 3 5 3 (k — k'). As mentioned in the previous sec- 
tion, the order parameter of the expansion is Hj 3 ^At, where At ~ \j\fmM. Using Eqs. 
(3.36) and (3.37), it can be estimated as, 



Hf'At-VnThV^l^Q 1 (3-38) 



2n 



AB C , (3.39) 
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at the resonance, k/a ~ m. Therefore, the higher-order terms are subdominant for the cubic 
interactions as long as AB^ < O(10 5 ). 

The quadratic interaction Hamiltonian (2.24) also gives corrections to Eq. (3.12). In a 
similar way as above, the order parameter of the corrections is estimated to be, 

tffAt-.yf. (3.40) 

Hence, we can use the perturbation safely when the feature in the power spectrum is negligibly 
small, AV^jV^ ~ q^m/H < 1, as discussed in the previous section. When this condition is 
not satisfied, we should solve the full equation of motion to get the mode function. In that 
case, the mode function contains the negative energy mode and the bispectrum has additional 
resonant contributions [32] . We can also estimate the cutoff scale A cut discussed in the last 
part of §2.2 by applying a similar argument to the derivative couplings with At = 2ir/m. 

3.4 Other Interactions 

In the previous subsections, we have seen that the bispectrum can be affected much by the 
resonance through the couplings K n and K<i at the specific scale. However, as can be seen in 
Eq. (3.2), their effects are suppressed by the Planck scale. This is because cubic interactions 
of the inflaton field are not contained in the couplings K n and K^- Then, they are induced 
through the gravitational couplings. Direct cubic interactions appear if there is the self 
interaction of the inflaton field, 

(3.41) 



4A^ 



They are also induced by the interactions that are higher order in A^, 

l ' : (9 X -^) 2 (^) 2 , (3.42) 



or violate the parity symmetry 



8A^ 



{dx-d^f. (3.43) 



Here, we estimate the magnitude of their effects on the bispectrum. The effects of the 
interactions (3.41), (3.42), and (3.43) on the background evolution and the power spectrum 
can be analyzed in a similar way as the couplings K&\ The interactions (3.41) and (3.42) 
have a little effect because they are higher order in the g-parameters, while the interaction 
(3.43) could have an effect as large as the couplings K<i with the g-parameter of the order of 
Q( € v f\t=o)^ ■ I n the parity- vaiolating case, the kinetic mixing term (dx • dcp) could appear. 
We assume that this term is so suppressed that the background evolution and the power 
spectrum are not affected much. 

The interaction (3.41) induce the interaction Hamiltonian as, 

° 3 {¥) I d " x ^ dip? - (3 - 44) 
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Oscillatory components are induced in this Hamiltonian at the higher order in the q-parameters 
(see Eqs. (2.10) and (2.11)). 

On the other hand, the interactions (3.42) and (3.43) induce the interaction Hamiltonian 



as. 



*hX<t> 



d 3 x ip 



2ip 2 



r(V^) 5 



(3.45) 



and 



4A^ 



d 3 x tp(d<pf 



(3.46) 



respectively. Using the dimensional argument in the previous subsection, their effects on the 
bispectrum can be estimated to be, 



AB 



(s) 



AB 



(h) 



q 2 (tvf x } =0 ) 



T7 ~ (ey/. 



x,t=o) 



AV, 



ABf ~ q{eyf-} = ^ 



X 



(evl 



x,t=o 



P< 
2 AV ( 



P, 



(3.47) 
(3.48) 



at the peak scale. Therefore, the corrections to the bispectrum through the interactions 
(3.41), (3.42), and (3.43) are enhanced respectively by the factors (m/H) 3 / 2 and (m/H) 2 , 
which can be much larger than the factor m/H for the couplings K n and K^. Since all mode 
functions appear with a derivative, these contributions scale as k^ 1 for finitely squeezed 
configurations where the subhorizon approximation (2.35) is appropriate. In the squeezed 
limit, on the other hand, we cannot use the subhorizon approximation and the bispectrum 
scales as k^ 3 as in the usual cases. 

If we require the conditions AP^ jP^ < 1 and m < (27r) 1 / 2 Arf for the perturbativity, 
they should be bounded from above as, 



AB ( C S \ ABf < O(10 5 ) x {evl 



Ah) 



X 



X 



AB 



(p) 



x,t=o) 

< O(10 6 ) x (ev/-^)' 



j_ 

' 14 



(3.49) 
(3.50) 



Therefore, the interactions (3.42) and (3.43) could induce large features in the bispectrum 
even when the features in the power spectrum are too small to be detected. Note that the 
scale of the resonance induced by the interaction (3.43) is different from that for the couplings 
Kd] The resonance occurs at K/a ~ m for the interaction (3.43) while K/a ~ 2m for the 
couplings Kj. If we consider much higher-order interactions, the scales of the resonance 
appear at the integral multiples of the mass scale. Though they are suppressed by the 
(/-parameters, they could still induce non-negligible features in the bispectrum. 



4 Summary 

Non-Gaussianity could contain various information on the physics behind inflation. In this 
paper, we have discussed the possibility that we could obtain hints on the heavy physics dur- 
ing inflation by analyzing local features in the primordial bispectrum. A heavy scalar field 
can leave non-negligible signatures in the primordial spectra through the parametric reso- 
nance between its background oscillation and the fluctuations in the inflaton field. We have 
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estimated the contributions from the resonance perturbatively by picking up the interactions 
relevant to it. The bispectrum is amplified at specific configurations, and its amplitude can 
be 0(1O 1-2 ), or as large as 0(1O 5 ) depending on the type of interactions within the param- 
eter region where the perturbative expansion is applicable. In particular, the resonance can 
give large contributions in finitely squeezed configurations, while the bispectrum cannot be 
large in the squeezed limit unless the modification to the power spectrum is large as in the 
case that the consistency relation is satisfied. Since we cannot observe the exact squeezed 
limit, these contributions could practically be most important in actual observations of the 
bispectrum in that limit. We have also found that there is a relation between the scales at 
which the features appear in the bispectrum and the power spectrum, and that the feature 
in the bispectrum can be much larger than that in the power spectrum. Though we have not 
discussed the observability of localized features in the bispectrum, it will be discussed in the 
upcoming paper by one of the present author [46]. If the correlated features are observed 
at the characteristic scales in the primordial spectra, it will indicate the presence of heavy 
degrees of freedom. By analyzing them, we could obtain some information on the physics 
behind inflation. 
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A Interactions Relevant to the Resonance 

Here, we discuss which cubic interactions are most relevant to the resonance. The cubic 
interactions for the action (2.1) have been obtained in Refs. [47, 48]. They can be separated 
into the gravitational and matter parts as, 



= -a 3 MpH 2 J d 3 x 3a 2 + ^«V 2 /3 + (V 2 /3V 2 /3 - ViV^VV/3) 
H$ = -a 3 J d 3 x P u (Xl J + al 2 ") + j, K L (2X{ J + aX{ J )X« L 

+ P u , K {Xl J + aX{ J )ip K + ^P/.jc^V + \pij,kl,mnX{ J Xf L X™ N 

2 D 

+ -^PlJ,KL,M X l X l <P + 7;PlJ,K,LX 1 <p <p + -Fj jK<p If if 
I 2 6 



a, (A.l) 



(A.2) 



Here, a and (3 are respectively the perturbations in the lapse and shift functions at the first 
order, which are given by solving the constraint equations as, 



V 2 /? 1 
-3a + 



a 2 H 2M^H 2 



PiV 1 - Pij,k^V<P K + (Pi j + Pu,LM^ L ^ M )^ L 4> M a - </>V) • 

(A.4) 
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The function X^ J (n = 1,2,3) represents the perturbations in X IJ = — {d ll (j) I ){d^(t) J )/2 at 
the n-th order. They can be read from 



^lE^-iv^W, (A.5) 



2(1 + a) 2 ^ n 2a 2 

^ ' 71=0 



where 



V IJ = <^ J , (A.6) 
Vl J = 2<j>V<p J \ (A.7) 

Vj J = ip V J - 4 V W (/ V iV > J) , (A.8) 

V 3 IJ = -\v i p^ I Vw J \ (A.9) 
cr 

v 4 IJ = ^(v^vVxv^vv 7 )- (A. 10) 

Here, () in the superscript represents the symmetrization. 

Since the resonance occurs deep in the horizon, the interactions with more derivatives 
are more relevant to the resonant enhancement of the bispectrum. Taking the terms with 
more derivatives, the lapse and shift functions are estimated to be, 

V 2 (3 [ev _}P_%± 2 ^ dl + g^)e~ 2Ft sm 2 (mt) 

a 2 H^y 2 M p H z 2 ' 1 ' 

at the leading order in the slow-roll parameter ey and the fraction f x . Though the energy 
fraction of the heavy scalar field f x contains an oscillatory component, 5 

/x^H fx,t=oe- 2Tt sin(2mf), (A. 13) 

we have neglected it because the decay rate T is assumed to be much smaller than the mass 
scale m, and then it cannot lead to a large enhancement. 

Counting the number of derivatives on (p, we can find the (<9(/?) 3 -type interactions in the 
terms 

qPij,kl,mnX[ j X^ L Xf IN , Pi.j,klXI j Xf L . (A. 14) 

However, these terms are absent unless we include the interaction (<90) 4 , or those that are 
higher order in or violate the parity symmetry (j) — > —(f). In §3.4, we discuss the resonant 
enhancement of the bispectrum induced by the (3</?) 3 -type interactions. 

The next candidates are provided by the (/?(3(/?) 2 -type interactions. They can be found 
in the terms 

_L_(V 2 /3V 2 /3- V.V^VV^a, P IJtK Xi J ^ K , (A.15) 
~P JJ>KL (2X 2 7J + aX( J )X« L , Pjj(Xi J + aX{ J ). (A.16) 



3 In the case that Y <C H, the suppression factor should be replaced by H/m. 
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The first term is higher order in the slow-roll parameter ey and the second term is absent 
unless parity-violating interactions like x L P{d L P) are introduced. Hence, the last two terms 
(A. 16) provide the cubic interactions that are most relevant to the resonance when only the 
couplings K n and Kd exist. In terms of zs and c s , these terms can be explicitly written as, 



^relevant ~ Z ^ I ^ X 



r 2 



(3 -2c^ + -|(V^ 



1 



+ — V^ViW } (A. 17) 



M p V 2 ./ | 2 



„2 



(3 - 2 C 2 )^ + %V^) 2 



a- 



+ [z\ + 2( %1 + g d2 )e~ 2ri sin 2 (mt)] VXV^ViW j . (A.18) 
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